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I. INTRODUCTION 



To explain the absence of the CP symmetry violation in interquark interactions Peccei and 
Quinn |l| suggested that a new symmetry must be present. The breakdown of this gives rise 
to the axion field proposed later by Weinberg and Wilczek j^]. And it was Wilczek who 
presented the first analysis of possible effects caused by axions in electrodynamics [4J. Notice 
;hat the idea to include an extra pseudoscalar field into electrodynamics was proposed by Ni 
5! as early as 1974, and so this date can be treated as the birth year of the axion prototype. 

Axions belong to the main candidates to form the dark matter, see, e.g. js] and references 
cited therein. New arguments for the materiality of axion theories were created in solid states 
physics. Namely, it was found recently |7J that the axionic-type interaction terms appear in 
the theoretical description of a class of crystalline solids called topological insulators. Axion 
electrodynamics gains plausibility by results of Heht et al [8[ who extract the existence of a 
pseudoscalar field from the experimental data concerning electric field-induced magnetization 
on Cr 2 03 crystals or the magnetic field- induced polarization. In other words, although their 
existence is still not confirmed experimentally axions are stipulated at least in the three funda- 
mental fields: QCD, cosmology and condensed matter physics. 

There are many other interesting aspects of axion electrodynamics. In particular, its reduced 
version (corresponding to the external axion field linear in independent variables) was used by 
by Carroll, Field and Jackiw (CFJ) [9] to examine the possibility of Lorentz and CPT violations 
in Maxwell's electrodynamics. In addition, just the interaction Lagrangian of axion electrody- 
namics generalizes the Chern-Simons form e abc A a V b A c [uj to the case of (l+3)-dimensional 
Minkowski space. 

Let us present more arguments for materiality of axion electrodynamics which are very 
inspiring for us. Recently new exactly solvable models for neutral Dirac fermions had been 
discovered 11]. j^j ]. These models involve the external electromagnetic fields which do not solve 
Maxwell equations with physically reasonable currents. However, these fields solve equations 
of axion electrodynamics. We had classified exactly solvable quantum mechanical models with 
matrix potentials [l3j|, [l^, and superintegrable models of cold neutrons 15]. Some of these 
systems also include external fields which solve equations of axion electrodynamics. In addition, 
these field equations appear to be a relativistic counterpart of Galilei invariant systems classified 
in [lf|. Thus we have a particular interest to study equations of axion electrodynamics, and 
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theory. However, except the analysis 
~y\ , we do not know any systematical 



we will do it using the tools of group theory. 

Group theory, and especially the theory of Lie groups is one of the corner stones of modern 
theoretical physics. Symmetries of Lagrangians and of the corresponding motion equations 
form a very essential constituent part of any physical 
of symmetries of the CFJ model presented in paper 
investigation of symmetries of axion theories. Notice that such an investigation would generate 
group-theoretical backgrounds for axion models and enable to construct their exact solutions. 

In the present paper we make the group classification of the field equations of axion elec- 
trodynamics with arbitrary self interaction of axion field. The considered model includes the 
standard axion electrodynamics as a particular case. We prove that an extension of the basic 
Poincare invariance appears only for the exponential, constant and trivial interaction terms. 
These and other results of group classification are presented in Section 3 and Appendix A. 

In addition, we carry out the group analysis of two other theories which are close to axion 
electrodynamics. Namely, we describe Lie symmetries of the field equations of classical electro- 
dynamics modified by adding the Chern-Simons term, and symmetries of the CFJ model. As 
it is shown in Appendix B, the maximal continuous group of Chern-Simon electrodynamics is 
the 17-parametrical extended conformal group. 

A special subject of our analysis are conservation laws which correspond to found symmetries. 
They are discussed in Section 4, where we present a simple proof that the interaction between 
the electromagnetic and axion fields does not affect the energy-momentum tensor. 

In Section 5 we present selected invariant solutions of field equations of axion electrodynam- 
ics. Some of these solutions play the key ro 
quantum mechanics in both relativistic llf . 



e in formulation of exa ctly solvable problems of 



12[ and nonrelativistic [15(, [18] approaches. 



In Section 6 we analyze plane wave solutions which are smooth and bounded functions 
which generate positive definite and bounded energy density. We show that these solutions 
describe waves whose group velocity can be superluminal. Nevertheless, they are causal since 
the corresponding energy velocities are smaller than the velocity of light. 

An important constituent of any relativistic model is its nonrelativistic limit. This is true also 
for models including massless fields. As it was shown long time ago 19| , there exist a reasonable 
(and very important) nonrelativistic approximation for the Maxwell equations, which makes 
them invariant w.r.t. the Galilei group. Namely, in this approximation we obtain equations 
of Faraday electrodynamics. This result justifies Galilei invariance of quantum mechanical 
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systems including particles interacting with an external electromagnetic field. 

A natural question arises whether it is possible to extend this result to the case of field 
equations of axion electrodynamics. Notice that the correct definition of the nonrelativistic 
limit of a physical model is by no means a simple problem in general and in the case of theories 
of massless fields in particular, see, for example, [20[. Such limit is not necessary unique, and 
simple passing the speed of light to infinity we can obtain a physically meaningless theory. 

By definition, any relativistic system is invariant w.r.t. the Poincare group P(l,3), and a 
correct nonrelativistic approximation of this system should be invariant w.r.t. the Galilei group 
G(l,3). Thus to obtain a well defined nonrelativistic limit it is necessary to take a care on the 



attending trans: 
and Wigner 



'ormation P(l,3)— >-G(l,3). This idea had been proposed long time ago by Inonii 
2l| who presented definitions and justifications for such transformation. It is a 
special limiting procedure called contraction, which is an important subject of modern group 
theory. 

In Section 7 we find a nonrelativistic limit of equations of axion electrodynamics with using 
a generalized Inonii- Wigner (IW) contraction. As a result we prove that the Galilei-invariant 



16|, are nothing but a 



wave equations for an abstract ten-component vector field, deduced in 
contracted version of the field equations of axion electrodynamics. 

Appendix A includes a rather detailed proof of the results formulated in Section 3. Fi- 
nally, in Appendix B we present the results of group analysis of the field equations of classical 
electrodynamics modified by adding the Chern-Simons terms, and of the CFJ model. 



II. FIELD EQUATIONS OF AXION ELECTRODYNAMICS 

Let us start with the following model Lagrangian: 

L = - \f»„F^ + j9F^ - V{9). (1) 

Here F^ is the strength tensor of electromagnetic field, F^ u = \e^ vpa F pa ', = 8^9, 9 is 
the pseudoscalar axion field, V(9) is a function of 9, k is a dimensionless constant, and the 

summation is imposed over the repeating indices over the values 0, 1, 2, 3. Moreover, the 
strength tensor can be expressed via four-potential A = (A , A 1 , A 2 , A 3 ) as: 

F V» = QV A v _ QV A H ( 2 ) 
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Setting in (PQ) — we obtain the Lagrangian for Maxwell field. Moreover, if 9 is a constant 
then (CQ) coincides with the Maxwell Lagrangian up to constant and four-divergence terms. Fi- 
nally, the choice V{6) = ^m 2 9 2 reduces L to the standard Lagrangian of axion electrodynamics. 

We will investigate symmetries of the generalized Lagrangian ([1]) with arbitrary V{9). More 
exactly, we will make the group classification of the corresponding Euler-Lagrange equations: 

d v F^ = K p v F"", (3) 
d v d v 9 = -|f^F> + F (4) 

where F = — In addition, in accordance with its definition, F^ v satisfies the Bianchi identity 

d v F» v = 0. (5) 

Substituting fl2]) into (j3J) one obtains the second order equation for potential A^. 

d v d v A» = -k Pu F^ (6) 

provided satisfies the Lorentz gauge condition: 

dpA» = 0. (7) 

Just the system of equations (EJ) ([5j) will be the main subject of group classification. In 
addition, we shall discuss symmetries of field equations of the Chern-Simons electrodynamics, 
i.e., of the system including equations (j3J) and (jSJ). In this theory is treated as an external 
field whose motion equation is not specified, i.e., equation @ is ignored. 



III. GROUP CLASSIFICATION OF SYSTEMS © © 

Equations Q-Q include an arbitrary function F(6) so we can expect that the variety of 
symmetries of this system depends on the explicit form of F. The group classification of these 
equations presupposes finding their symmetry groups for arbitrary F. 

In this section we present the results of group classification while the related calculations 
details are given in Appendix A. 

The maximal continuous symmetry of system (jSjHdHJ) with arbitrary function F(8) is given 
by Poincare group P(l,3). On the set of solutions of equations ^3J) — (JHJ) written as a seven 
component vector 

^ = column(F°\ F 02 , F 03 , F 23 , F 31 , F 12 , 9) (8) 



infinitesimal generators of this group take the following form: 



PfJ, C^l) JflU %fldv Xyd^i ~\- 

where indices fi and v independently take values 0, 1, 2, 3, 



(9) 



S 



ab 



£-abc 



\ 



\ 







/ 



.5, 



Oc 



Sr 



\ 







a, b, c 7^ 0, 



S c are 3x3 matrices of spin 1, whose entries are (S c ) a b = ie ca b, e ca b is the Levi-Chivita symbol, 
and the dots denote the zero matrices of an appropriate dimension. 

Alternatively, spin matrices S^ can be represented as the first order differential operators: 

S^v = F^\dp\v — F v \d F \^. (10) 

Such notation is both convenient and usual for group analysis of differential equations. 

Operators (Q form a basis of the Lie algebra p(l,3) of the Poincare group. 

For some special functions F(9) symmetry of system ([3])-([5]) appears to be more extended. 
Namely, if F = 0, F = c or F = bexp(a8) then the basis ([§]) of symmetry algebra of this system 
is extended by the following additional operators P 4 , D and X: 

p 4 = d e , D = x d + xA - \F" v d F ^ if F(6) = 0, 

PA = d e if F(9) = c, (11) 

X = aD - 2P 4 if F(9) = be a0 . 

Operator P4 generates shifts of dependent variable 9, D is the dilatation operator generating a 
consistent scaling of dependent and independent variables, and X generates the simultaneous 
shift and scaling. Note that arbitrary parameters a, b and c can be reduced to the fixed values 
a = ±1, b = ±1 and c = ±1 by scaling dependent and independent variables. 

Thus the continues symmetries of system (J3j) — (JHI) where F{9) is an arbitrary function of 9 are 
exhausted by the Poincare group. The same symmetry is accepted by the standard equations 
of axion electrodynamics which correspond to F(9) = —m 2 9. In the cases indicated in ( 1TT1) we 
have the extended Poincare groups. 

Notice that symmetries of the equations §6§, ([7]), ([5]) for potentials also are generated by 
infinitesimal operators of the form ([9]) where 

S ltt , = A lt d Av -A u d All . (12) 

Additional symmetries again are given by equations (jlip where, however, D — > x M <9 M — A^8a^- 



IV. CONSERVATION LAWS 



An immediate consequence of symmetries presented above is the existence of conservation 
laws. Indeed, the system flHJ-flSJ) admits a Lagrangian formulation. Thus, in accordance with 
the Noether theorem, symmetries of equations (J3j) ([5]) which keep the shape of Lagrangian ([!]) 
up to four divergence terms should generate conservation laws. Let us present them explicitly. 

First we represent generators , (1121) and (TTT|) written in terms of the variational variables 

and A 4 = 9 in the following unified form: 

Q = ^d, + ^d A -r (13) 

where the summation is imposed over the values r = 0, 1, 2, 3, 4 and \x = 0, 1, 2, 3 . 



Conserved current corresponding to symmetry (I13p can be represented as 221 ]: 



J -=^m^ +rL - ra ^m^y (14) 

The basic conserved quantity is the energy-momentum tensor which corresponds to symme- 
tries P M presented in (jUj). In this case 

ip T = and ^ = 1 (15) 

where \x successively takes the values 0, 1, 2, 3 . Substituting ([Tj), ffl5]) into (fT4"]l and using three 
dimensional notations 



Oa 



E a , F ab = e abc B c . (16) 



we find the components of the conserved energy-momentum tensor in the following form: 

T 00 = i(E 2 + B 2 + p 2 + p 2 ) + V(9), T 0a = T a0 = e abc E b B c + p°p a , 

2 j (17) 

T ab = _ E a E b _ B a B b + p a p b + gab ^.2 + g2 + ^2 _ p 2 _ 2 y(0)). 

The tensor T M!/ is symmetric and satisfies the continuity equation dJT^ = 0. Its components 
T 00 and T 0a are associated with the energy and momentum densities. 

It is important to note that the energy-momentum tensor does not depend on parameter k 
and so is not affected by the term jdF^F^" present in Lagrangian ([1]). In fact this tensor is 
nothing but a sum of energy momenta tensors for the free electromagnetic field and scalar field. 
Moreover, the interaction of these fields between themselves is not represented in (IT7|) . 



The conservation of tensor (|T7j) is caused by the symmetry of Lagrangian ([T]) w.r.t. shifts of 
independent variables x M . The symmetries w.r.t. rotations and Lorentz transformations give 
rise to conserving of the following tensor: 

G av " = x a T» v - x v T^ (18) 

which satisfies the continuity equation w.r.t. the index fj,. In particular, for a, v = 1,2,3 
equation ( fTBl) represents the angular momentum tensor. 

The tensors (TTTI) (HTH]) exhaust the conserved quantities whose existence is caused by the Lie 
symmetries of equations with arbitrary function F(9). 

The additional symmetries presented in (11 II) are neither variational nor divergent symmetries 
of Lagrangian ([1]), and so they do not generate conservation laws. However, we can indicate 
another conservation laws which have nothing to do with Lie symmetries. 

First let us note that equation ((3]) in itself can be written in the divergence form d v Tv = 0, 
where 

T^ = - kOF^ 

is the antisymmetric conserved tensor. In addition, this equation can be rewritten as = 
where 

f = d u F^ - k Vv W v (19) 

is a conserved current. Changing equation ([3]) by ( TT9l) with ^ we obtain the system which 
represents the field equations of axion electrodynamics with nontrivial currents. 
Equation (0J in its turn can be represented as 

= F{6) 

where 

J» = I f + kF^A u . (20) 

If F = then current (1201) satisfies the continuity equation. 

In addition to ( |T71) -(l20l) there exist the infinite number of (trivial) conserved currents cor- 
responding to the gauge symmetries of Lagrangian ([1]). An example of such conserved current 
is: 

ATM = FV Puip (e) (21) 

where <p{0) is an arbitrary differentiable function of 9. Vector satisfies the continuity 
equation d^N^ = provided equations ([3]) are satisfied (remember that p„ = d v 6). 
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V. SELECTED EXACT SOLUTIONS 



The field equations of axion electrodynamics form a rather complicated system of nonlinear 
partial differential equations. However, this system admits an extended symmetry algebra 
which makes it possible to find a number of exact solutions. Here we present some of these 



solutions while the completed list of them can be found in 



23]. 



The algorithm for construction of group solutions of partial differential equations goes back 



to Sophus Lie and is expounded in various monographs, see, e.g., [22|. Roughly speaking, to 
find such solutions we have to change the dependent and independent variables by invariants 
of the subgroups of our equations symmetry group. Solving equations (J3j) ([5]) it is reasonable 
to restrict ourself to three-parametrical subgroups of P(l,3) which enables to reduce (E])-© to 
systems of ordinary differential equations. The complete list of these subgroups can be found 
in [24] . 

To make solutions of equations (EI) ([5j) more physically transparent, we write them in terms 
of electric field E and magnetic field B whose components are expressed via the strengths tensor 
as shown in f|T6|) . In addition, we rescale the dependent variables such that k — > 1. 



A. Plane wave solutions 



Let us present solutions of system (J3j) (J5j) which are invariant w.r.t. subalgebras of p(l,3) 
whose basis elements have the following unified form: (P 1; P 2 , kPo + eP 3 ) where e and k are 
parameters satisfying e 2 ^ k 2 , while Pi, P 2 , P3 and Pq are generators given in (CQ). 

The invariants u of the corresponding three-parametrical group should solve the equations 

P 1 cu = 0, P 2 u = 0, (kP + eP 3 )uj = 0. (22) 

Solutions of ( 1221) include all dependent variables E a ,B a ,9 (a = 1, 2, 3) and the only independent 
variable u = exq — kx 3 . Thus we can search for solutions which are functions of u only. As a 
result we reduce equations (J!])-©, (TT6l) to the system of ordinary differential equations whose 
solutions are 

B>i = -kci9, B 2 = ec 1 + kc 2 , B 3 = c 3 , 

, , (23) 
Ei=ec 2 + kci, E 2 = ec 1 9, E 3 = c 3 9 - c 4 (e : - k ) 
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where C\,...,C4 are arbitrary real numbers. The corresponding bounded solutions of equation 
(TjJ with F = -m 2 9 are: 

9 = a M cos fiu + r M sin fiu + (24) 

where a M , r M and \i are arbitrary constants restricted by the following constraint: 

2 f 2 , cj + m 2 ^ 
/i= ^ + ^3pJ- (25) 

Notice that for the simplest nonlinear function F = X8 2 equation (jlj) is reduced to Weierstrass 
one and admits a nice soliton-like solution 

= t an h 2 (co + C) (26) 

where C is an integration constant. The related parameters e, k and A should satisfy the 
conditions 

c 2 

e 2 = k 2 + —2-5, Ac 3 c 4 = 12. (27) 

o 

The corresponding magnetic, electric and axion fields are localized waves moving along the 
third coordinate axis. 

In analogous (but as a rule much more complicated) way we can find solutions corresponding 
to the other three dimensional subalgebras of the Poincare algebra. One more and rather specific 
solution of equations Q-© with k — 1 and F = (obtained with using the subalgebra spanned 
on basis elements (J12 + kP + ePi, P2, P3)) can be written as follows: 



Ei = e(ck sin(u;) — dk cos(a;)), E 2 = e(ck cos(a>) + dk sin(oj)), E 3 = e, 

k k 
Bi = — E 2 , B 2 — -E\, B 3 = 0, 9 = axo + ux 3 + fi 

e e 



(28) 



where e, Ck, dk, s, k, a, is, \x are constants satisfying the following conditions: 

e 2 - k 2 = ve - ak, e^O. (29) 



Solutions ff28|) depend on two different plane wave variables, i.e., u = ex — kxi and axo+vxi. 
They satisfy the superposition principle since a sum of solutions with different e,k,Ck and dk 
is also a solution of equations (jSJ) - © with k — 1 and F — 0. Thus it is possible to generate 
much more general solutions by summing up functions (1281) over k and treating and dk as 
arbitrary functions of k. 
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B. Radial and planar solutions 



Consider solutions which include the Coulomb electric field. They can be obtained using 
invariants of the subalgebra spanned on (Ji 2 , J23, J31) arid have the following form: 

where ip is a function of xq and r = \fx\ + x\ + x\ satisfying the following equation: 

d 2 cp d 2 (p (c\ 2 \ cic 2 . . 

n + m ( 31 ) 



<9r 2 cte 2 , yr 4 / r 3 

Setting in ( )30l) ci = we come to the electric field of point charge which is well defined 
for r > 0. A particular solution for (|3T|) corresponding to C\ = —q 2 < and C2 = is 
97 = c 3 r sin(mxo)e~~ which gives rise to the following field components: 

Q 2 X(i q 2 9x^ <p 

B a = — — > = -3-^, = c 3 sm(mx )e~~. (32) 

The components of magnetic field B a are singular at r = while E a and 9 are bounded for 
< r < 00. 

Separating variables it is possible to find the general solution of equation (I3TI) . see [23]. 
One more solution of equations (j3J)-(jSJ) for F = with a radial electric field is: 

E a = X f 2 . (33) 

The corresponding magnetic and axion fields take the following forms: 

rj ^1^3 R X 2 X 3 R X / X 

By = B 2 = —5—, B 3 = — -, 9 = arctan — 

r x r x r A \x 3 



where x = \/x\ + x\. 

The electric field ( 133]) is requested in the superintegrable model with Fock symmetry pro- 
posed in [12]. 

Let us present planar solutions which depend on spatial variables x\ and x 2 . Namely, the 
functions 

E x = x x ( Cl x C3 - 2 + c 2 x- 2 - c A , E 2 = x 2 ( Cl x C3 " 2 + c 2 x~ 2 - C3 ) , E 3 = 0, 

(34) 

B x = x 2 ( Cl x C3 ~ 2 - c 2 x~ 2 - C3 ) , B 2 = Xl (c 2 x~ 2 - C3 - cix C3 ~ 2 ) , B 3 = 0, 

9 = c 3 arctan — + C4 (35) 

X\ 
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where ci, c 4 are arbitrary parameters, solve equations d2|)-(J5]) with k = 1 and F = 0. 
In particular, for C2 = 0, c 3 = 1 and C2 = 0, C3 = — 1 we have: 

E 1 = -B 2 = C ^ B 1 = E 2 = ^, B 3 = E 3 = 0, 6 = arctan^ (36) 

and 

^ = -^ = ^1 5l = £ 2 = ^, S3 = i?3 = 0) ^ = arctan ^. (37) 

Solutions (IMjl . (1351) can be found with using invariants of a subgroup of the extended Poincare 
group whose Lie algebra is spanned on the basis {Po, P 3 , Ju + P4), see equations (TTTj) for 
definitions. 

Let us present an example of solutions describing fields in a constantly charged space. The 
related equation (|3]) for /i = should be changed to 

VE = pB+j 

(compare with ( 1461) ) where jo represents a constant charge density. The remaining equations 
(EI)— (EJ) are kept uncharged, and the considered system is solved by the following fields: 

Ei = - j xi In x, E 2 = In x, E 3 = B 3 = 0, 

Bi = ^j x 2 Unx- B 2 = -^j xi Hnx - ^ , (38) 

x 2 

= 2 arctan — . 

Xi 

The fields fl36l) and fl38l) appears in superintegrable models for particles with spin 1 and | 
correspondingly [18[ while the fields (1371) are requested in the exactly solvable system described 



by the Dirac equation 111 ]. 



VI. PHASE, GROUP AND ENERGY VELOCITIES 

In this section we consider some of the found solutions in more detail and discuss the 
propagation velocities of the corresponding fields. There are various notions of field velocities, 
see, e.g., 25] j^J, 2^]. We shall discuss the phase, group and energy velocities. 

Let us start with the plane wave solutions given by equations (123]) and ( |24|) . They describe 
oscillating waves moving along the third coordinate axis. Setting for simplicity c 2 = c 3 = C4 = 
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V* = -EH = ~ ( 41 ) 



= we obtain: 

B 1 = Cl k9, B 2 = -cie, B 3 = 0, 

(39) 

Ei = —cik, Ei = —c\eQ, E% = 0, 6 = a M cos(//(££o — A;x 3 )). 

Here e, k, and a M are arbitrary parameters which, in accordance with (|25l) . should satisfy the 
the following dispersion relations: 

{e 2 - k 2 ){^ 2 - c\) = m 2 . (40) 

If m 7^ the version /i 2 = c\ is forbidden, and we have two qualitatively different possibilities: 
fi 2 > c\ and yU 2 < c\. 

Let /i 2 > c 2 then (e 2 — k 2 ) = jjrz^ > 0. The corresponding group velocity V g is equal to the 
derivation of s w.r.t. k, i.e., 

de k 
dk e 

Since e > k, the group velocity appears to be less than the velocity of light (remember that we 
use the Heaviside units in which the velocity of light is equal to 1). 

On the other hand the phase velocity V p — | is larger than the velocity of light, but this 
situation is rather typical in relativistic field theories. 

In the case /i 2 < c 2 the wave number k is larger than e. As a result the group velocity (|4ip 

exceeds the velocity of light, and we have a phenomenon of superluminal motion. To understand 

wether the considered solutions are causal let us calculate the energy velocity which is equal to 

the momentum density divided by the energy density: 

7->03 

V * = too- (42) 
Substituting ( 139]) into (fT7|) we find the following expressions for T 00 and T 03 : 

T 00 = £ (£ 2 + k 2^ + i m 2^ T 03 = £k ^ 

where $ = c 2 (6 2 + 1) + /i 2 (a 2 - 6 2 ). Thus 

e ~ {e 2 + P)$ + |m 2 ^ 2 < e 2 + k 2 < ' 
and this relation is valid for e > k and for e < k as well. 

We see that the energy velocity is less than the velocity of light. Thus solutions ( 13"9"j) can be 
treated as causal in spite of the fact that for /z 2 < c 2 the group velocity is superluminal. 
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Analogously, analyzing dispersion relations (T2"9"1) we conclude, that the electromagnetic fields 
291) propagate along the third coordinate axis with the group velocity 



V 9 



dk 



' (43) 



VTT5 



were 



u 2 — a 2 



(44) 



(2k -a) 2 ' 

If — 1 < 5 < the group velocity (144 j) is lager than the velocity of light. However, the 
corresponding energy velocity (142 p which is equal to 

2e/c + 2z/a 
+ k 2 + v l + a" 1 

cannot exceed the velocity of light since 2ek < e 2 + k 2 and 2ua < v 2 + a 2 . 

Thus solutions (128]) also can propagate with a superluminal group velocity. However, these 
solutions are causal since their energy velocities ( 1431) are subluminal. Analogous results can be 
proven for soliton-like solutions ( l26l) . 



VII. NONRELATIVISTIC LIMIT 

To find a nonrelativistic limit of the field equations of axion electrodynamics we shall use 
the generalized IW contraction which guaranties Galilean symmetry of the limiting theory. 

Let us consider more general field equations including currents (see ( DSD ), but restrict our- 
selves to the limiting case of the zero axionic mass. More explicitly, we start with the following 
system: 

V-E = «p-B+j , (46) 
<9 E- V x B = n(p B + p x E) +j, (47) 
d B + V x E = 0, (48) 

V ■ B = 0, (49) 
d oPo - V • p = -«E • B + j 4) (50) 
<9oP-Vp = 0, (51) 

V x p = 0. (52) 
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If Jo — J4 — and j = then the subsystem (SSD (SHJ) reduces to equations (j3J), ©, while 
the subsystem (I50l)-(l52l) becomes equivalent to (jnj) (remember that po = ^o^ 1 and p = V0). 

Like ([3])-([5]) the system with currents, i.e., ( H6|) - (l52l) is Poincare invariant. The related 
representation of the Lie algebra of Poincare group can be obtained by the prolongation of the 
basis elements ([9]) to the first derivatives of 9 and adding analogous terms acting on components 
of the current four-vector j = (j ,j) = (j ,ji,j 2 ,h)' 

Po — do, P a = d a , 

•L = x a d b - x b d a + B a d Bb - B b d B a + E a d Eb - E b d E a + P a d pb - P b d r + fd jb - j b d ja , (53) 

Joa = x d a + x a d + Babe (E b d B c - B b d E o) + p°d pa - p a d p0 + j°d ja - j a d j0 . 

Generators (1531) do not include differentials w.r.t. j 4 since this current component is not changed 
under Lorentz transformations (i.e., it should be scalar). 

Being applied to basis elements of algebra p(l,3) the IW contraction consists of the trans- 
formation to a new basis 

Jab J'ab = Jab, Joa ~* J'o a = £ Joa, Po — * Pq ~ £ 1 pQi Pa ~^ P' a = Pa (54) 

where e is a small parameter equal to the inverse speed of light. In fact we deal with the gen- 
eralized IW contraction since Pq is proportional to the inverse power of the small parameter. 
In addition, the dependent and independent variables in fl53|) undergo the invertible transfor- 
mations E a — > E' a ,B a — y B' a ,p fl — y p'^, x M — y x'^ where the primed quantities are functions 
of all the unprimed ones and of e. Moreover, the transformed quantities should depend on the 
contracting parameter e in a tricky way, such that all transformed generators P' J' ab and J' Qa 
are kept nontrivial and nonsingular when e — y 0. 

The contractions of relativistic bi-vector fields (like E, B) and four- vectors p^ has been 



described in papers [28j and (29|. However, we need to contract simultaneously two subjects, i.e., 
the basis elements of the Lorentz algebra given by equations fl53|) and the system of equations 
(|46l) - (l52|) . which is a much more sophisticated problem. 

In order that transformation ( 1541) be nonsingular in e it should to be attended by the following 
transformations of the dependent and independent variables: 

Xq t £2^0) -^a, 

e 
2' 



Po=Po, p' = §(E + p), E^e-^E-p), B' = B, (55) 



fo = e 1 Uo + k), n = |(io - 34), j' = j. 
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Making changes (|55[) in equations (|53p . applying transformation (|54p and tending e — > 
we obtain the following set of first order differential operators which for a basis of the Galilei 
algebra: 

Pq — dt, P' a = d a , 

J' ab = x a 8 b - x b d a + B' a d B , b - B' b d B , a + E' a d E , b - E' h d E , a 



+ p' a d p , b - p' b d p/a + j' a d fb - j' b d f a, 



(56) 



J' Qa = td a + e abc (p ,b d B , - B /b d E ,c) + p'°d E , a - p' a d pl0 - j /a d f o. 
To obtain a consistent system of equations, we have to make changes (155]) not directly in 
equations (|46|) - (l52|) . but in the equivalent system which includes equation (1471) . (f49|) . (!51~|) and 
sums and half divergences of pairs of equation (jlBT) . ( l50i) and (HE]), ( |52i) . Then equating terms 
with lowest powers of e we obtain the following system: 

d t p' + V-E' + kB>-E' = f Q , 
d t p' - V x B' - k( P ' B' + p' x E') = j, 
V • p' — up' ■ B' = j' 
V-B' = 0, 
d t B' + V x E' = 0, 
d t p' -Vp' = 0, V x p' = 

and p'° = d t 9', p' = V6'. 

Just equations (l57j) present the nonrelativistic limit of system (j46l) - (!52l) . The Galilei invari- 
ance of system (1571) can be proven directly using the following transformation laws which can 
be found by integrating the Lie equations for generators (1561) : 

X — 7- x + vt, t — > t, 

Po^Po-v-p', p'^p\ B' = B'+vxp', (58) 
E'4E'-vxB' + + v(v ■ p) - ^v 2 p'. 
Thus we find the nonrelativistic limit for field equations of axion electrodynamics with zero 
axion mass and nontrivial currents. It is interesting to note that system ( l57j) has been discovered 



m paper 



16| , starting with the requirement of Galilei invariance for abstract vector fields. Thus 



our contraction procedure presents a physical interpretation for the Galilei invariant system for 



the indecomposable ten component field deduced in [16(, see equation (67) there. Namely, this 
Galilei invariant system is nothing but a nonrelativistic limit of the field equations of axion 
electrodynamics. 
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VIII. DISCUSSION 



The aim of the present paper is multifold. First we make group classification of field equations 
of axion electrodynamics which include an arbitrary function F depending on 6, and 

find the conservation laws generated by these equations. Secondly, we find all continuous 
symmetries of Chern-Simon electrodynamics and CFJ model which are closely related to axion 
electrodynamics. At the third place, we discuss exact solutions for axion and e.m. fields which 
can be obtained using found symmetries. Finally, we define a correct nonrelativistic limit of 
the field equations of axion electrodynamics. 

As a result of the group classification we prove that the Poincare invariance is the maximal 
symmetry of the standard axion electrodynamics and indicate the special forms of F for which 
the theory admits more extended symmetries. In accordance with ( ITT]) such an extension 
appears for trivial, constant and exponential interaction terms. 

In spite of the transparent Poincare invariance of the analyzed equations, these results are 
nontrivial, since possible extensions (jTT|) of the basic symmetries and the absence of other ones 
were not evident a priori. Our results form group-theoretical grounds for constructing of various 
axionic models. The detailed calculations of Lie symmetries of equations (jS))-© are presented 
in Appendix A. 

In addition, in Appendix B we describe Lie symmetries of classical electrodynamics modified 
by Chern-Simons term (i.e., symmetries of the subsystem ([3]), (jSJ) where p M is treated as an 
external field), and symmetries of CFJ electrodynamics. It is proven that the maximal contin- 
uous symmetry of the Chern-Simons electrodynamics is given by the 17-parametrical extended 
conformal group. 

Conservation laws and the corresponding symmetries of the CFJ model have been already 

studied in paper However, the symmetries which generate conservation laws are not the 

half of all Lie symmetries 22]. Our research justifies and completes the results of \y\ - In 

particular we prove the completeness of the list of space-time symmetries theory presented in 
I I 

paper [17j . and add this list by symmetry Di given by equation fIBlOl) . 

Analyzing conservation laws caused by equations (J3j) — ([5]) we prove that the energy- 
momentum tensor does not depend on interaction between the electromagnetic and axion fields 
and indicate that the additional symmetries ( TTTT) do not generate conservation laws. The first 
statement generalizes the observation present in book 30j to the case of arbitrary self interaction 
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of axionic field. 

Exact solutions of the field equations of axion electrodynamics are discussed in Sections 
5 and 6. There is a sufficiently large amount of group solutions for these equations, but we 
restrict ourselves only to those which can be potentially important to physical applications. 

Analyzing the plane wave solutions we recognize the possibility of the faster than light group 
velocities in axion electrodynamics. However the corresponding energy velocities are subluminal 
and do not lead to causality violation. 

Let us note that functions ( 1281) solve the field equations of CFJ electrodynamics which 
coincide with the system fl3]), fl5]) where p^ are constants. The existence of faster than light 
solutions for these equations was indicated in paper [9], the enhanced discussion of the plane 



31] . However, in CFJ electrodynamics the energy is 



wave solutions for ([3]), ([5]) can be found in 
not positive definite [9] while the energy density of axion electrodynamics, discussed in section 
5, is positive. In addition, in the CFJ theory variables p^ represent an external field while in 
our approach they are dynamical variables satisfying equations (Hj). 

An important aspect of the presented exact solutions is that they are requested in some 
problems o 



nonrelativistic and relativistic quantum mechanics. In particular, as it was indi- 
cated in [ill ], just the vectors of the electric and magnetic fields described by relations ( 1371) give 
rise to exactly solvable Dirac equation for a charged particle anomalously interacting with these 
fields. These fields are involved also into the superintegrable nonrelativistic system discussed 
in [15] . The electric field given by equation ( )33|) was used in [12| to construct a new exactly 



solvable QM model with Fock symmetry. The fields f )36|) and f )38|) are requested in superinte- 
grable models for particles with higher spins 18] . Using the formalism presented in 32] it is 
possible to construct also relativistic versions of these models. 

Summarizing, the presented exact solutions are not nice mathematical toys only. In contrary, 
they have a nontrivial physical content. These solutions give rise to exactly solvable QM systems 
enumerated in the previous paragraph, which in principle can be used to predict physically 
verifiable effects. That opens ways for finding new arguments for the real existence of the 
pseudoscalar axion field. 

Solutions discussed in Section 4 represent only a part of exact solutions which can be obtained 
using three dimensional subalgcbras of algebra p(1.3). The complete list of them can be found 
in 23]. Notice that among these solutions there are rather general ones which include six 
arbitrary functions 23]. 
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A special goal of this paper was to present a correct nonrelativistic limit of equations of 
axion electrodynamics. To achieve this goal we use the generalized IW contraction of the 
corresponding representation of the Poincare group. As a result we prove that the limiting case 
of these equations is nothing but the Galilei-invariant system for the ten-component vector 



field obtained earlier in paper 16] . This result casts light on the physical content of the model 



discussed in 



16| . In addition, the contraction procedure can be used for transforming solutions 



discussed in Sections 4, 5 (and other solutions found in [23j) to solutions of system 



Some results of this paper had been announced in conferences, see Proceedings [33 ]. 
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Appendix A: LIE SYMMETRIES of equations ©-© 

The group analysis of differential equations is a nice field of mathematics whose fundamentals 
had been created by Sophus Lie as long ago as in the end of the eighteenth century. The 
foundations of the group analysis are expounded in many books, the most popular of which is 



monograph 22j . Nevertheless we will present a very short sketch of the Lie algorithm addressed 
to more physically oriented readers. 

The basic idea of the classical Lie algorithm is to treat a differential equation (or a system of 
equations) as a manifold in the multidimensional Euclidean space whose basis elements include 
the dependent and independent variables and also the derivatives present in the equation. In 
our case we have four independent variables x°, x 1 , x 2 , x 3 , seven dependent variables F^ u , 9 and 
38 differential variables Fff = dxF^ , 6>„ = d v 9 and 9^ v = d^d v 9. Equations ©-© define a 
hypersurface T in the Euclidean space K49 with variables x 1 *, F^ u , 9, F^ , 9 tU , i(U „. Moreover, 
J 7 is a smooth manifold in 

Local continuous symmetries of system (J3J) ([SD can be treated as continuous transformations 
in which keep this manifold invariant. To find these symmetries it is possible to imply the 
infinitesimal invariance criterium, i.e., consider transformations close to the identity one: 

x » x >» = x k + F^ -> F'^ = F"" + erf v , 9 -> & = 6 + eo (Al) 
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where e is a small transformation parameter and rj^ u , a are some functions of x M , F^ v and 
9, and ask for the form invariance of equations (j3J) - flH]) w.r.t. the change of variables (IA1I) . 
Transformations (1A1[) can be formally represented as 

f-+f = (l+eQ)f 

where / is a vector whose components are independent variables and dependent variables 
F^ v , 9, and Q is the infinitesimal operator: 

Q = ed,+ l -^ u d F ^+ad e . (A2) 



Starting with (lAlj) we can find transformations for the differential variables: 

*7 *T = *T + , 9,^9*^ = 9^ + 8^, 9^^9*^ = 9^ + 8(7^ (A3) 

where 

rtf = D^rT) - FJiy^), <r„ = £>» - tf AJ D,(£ A ), 
^ = A,(<7„) - 0,,aA,(£ A ) 

and D M = a M + \F™d F v« + 9^d e + 0,„„fy„. 

Using ( 1A3I) . the invariance condition for system (J3l) — (J5j) can be written in the following form: 

Q { 2)F\r=o = (A4) 
where Qi-a is the infinitesimal operator ( 1A2I) prolonged to the first and second derivatives: 

Q(2) = Q + \)7 Q F7 + + \^ud 9 ^ (A5) 

and T is the manifold defined by relations (J3j) — (J5I) - 

Acting by operator (IA5I) on differential forms ©-([5]) and equating coefficients for linearly 
independent functions F^, 9 and their derivatives we obtain the following system of determining 
equations for the coefficients rf v and a: 

&* = o, tf = o, C + C = ^C> ( A6 ) 

o F ^ = 0, a ee = 0, (A7) 
^ + (a e - ^) (F + \f, v F^) - \rf v F^ - aF e = 0, (A8) 
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< = F^a u , e^ aaV a J = 0, 2a 0tll = d v d v ^ 

rf v - \e^ v a ^T + F^ae + F™^ a - F^ a = 0, (A9) 

+ Vf^v = 0' Vf^k = Vf^ v ' = F^tj-Ih 

where the subscripts F u \ and 6* denote the derivatives with respect to the corresponding vari- 
ables: £p vX = gp^s:, r]g U = etc., 5 V is the Kronecker symbol and there are no sums over 
the repeating indices in the last line of equation ( 1A9I) . 

In accordance with equations (1A6|) functions £ M do not depend on F^ v and 9 and, moreover, 
they are Killing vectors in the space of independent variables: 

= 2x ti f v x u - f»x v x v + (f v x v + dx» + (A10) 

where f 1 , d, e M and = — c 1 ^ are arbitrary constants. 

It follows from (1A7I) that a = tpi9 + <p 2 , where tp\ and tp 2 are functions of x^. Substituting 



this expression into (IA8[) we obtain the following equation: 

ifxOFe + <p 2 F e + (2£° - <pi) (f + ^F M „F>) 



(All) 

+ KTf v % v - Qd a d a ^ - d a d a ip 2 - 2p»d^ = 0. 



Let the terms 9Fg, Fg, F, and 1 be linearly independent. Then it follows from ( lAllj) and 



that (pi = ip 2 = 0, f v = and ifF^ = 0. Thus, using (jSOJ) and (lATOl) we obtain: 

^ = c ^ + e^, r]^ = (f a F u ct -c ua F >l a , a = 0. (A12) 

Substituting (jA"l2l) into (jA2]) we receive a linear combination of infinitesimal operators (Q 
which form a basis of the Lie algebra of Poincare group P(l,3). Thus the group P(l,3) is the 
maximal continuous symmetry group of system (TJ|)-([3]) with arbitrary function F{9). 

The possible extensions of this symmetry which appears for some particular functions F are 
enumerated in equations (11). 



Appendix B: Symmetries of Chern-Simons and Carrol-Field-Jackiw models 

Let us discuss symmetries of two models which are closely related to the axion electrody- 
namics. The first is the classical electrodynamics modified by adding the Chern-Simons terms. 
It is based on field equation (jSJ) , (JSJ) , but does not include the dynamical equation for p^ which 
is treated as an external field. 



21 



Let us suppose that vector p^ can be presented as a four-gradient of some scalar function 
9. Then, repeating the procedure used above we again come to the determining equations 
( 1A6|) . ( 1A7j) . ( 1A9j) while equation (1A7|) would be absent. Solving these determining equations 
we obtain the 17-dimensional symmetry algebra whose basis elements are given by equations 

and by the following equations: 

D 1 = F" u d F „, D 2 = x^-D 1 , 

(Bl) 

= 2 XfM D 2 

p < = I < B2 > 



where are operators defined by equation (flOl) . 

Operators (J9j) and ( 1B1[) form a basis of the Lie algebra of conformal group, while operators 
P4 and D 2 belong to the central extension of this algebra. 

Thus the field equations of Chern-Simons electrodynamics are invariant w.r.t. the conformal 
group and admit two additional symmetries, i.e., scaling of the strength tensor of electromag- 
netic field and shifts of function 9. 

If vector field p^ is not supposed a priori to be a four-gradient, we come to the 16-dimensional 
symmetry algebra of equations (EJ) , (J5j) . The explicit form of its basis elements can be obtained 
from Q) and (1B1I) by the change S^ v — > S^ u and D 2 — > D 2 where 

Sfw = + p^dpv - p u d p „ , D 2 = D 2 - p^dpv . 

The second model which we consider here is the Carrol-Field- Jackiw one 9j. This is a 
particular version of the Chern-Simons electrodynamics corresponding to a constant vector p M . 
Considering equations ([3]), ([5]) with the additional condition 

d v p„ = 0, (B3) 

we conclude that in this case the symmetry algebra is reduced to p(l, 3) © Ai where p(l, 3) is the 
extended Poincare algebra whose basis elements are operators ([9]) and operator D 2 , while Ai 
is the one-dimensional algebra spanned on D\. The reason of this reduction is that equations 
(1B3j) are not invariant w.r.t the conformal transformations generated by operators K^. 

Thus the system ([3]), ([5]), (IB3I) is not invariant w.r.t. the conformal transformations, but it is 
still invariant w.r.t. the extended Poincare group P(l,3) which includes the ordinary Poincare 
group and accordant scalings of dependent and independent variables. Equations ( 1B3j) require 
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that p^ are (arbitrary) constants. If these constants are fixed, the symmetry group P(l,3) 
is reduced to its subgroup which keeps them invariant. Namely, for time-like, space-like and 
light-like vectors p = (po,Pi,P2,P3) we can set 

p= Ox, 0,0,0), (B4) 

P= (0,0,0,/i), (B5) 

and 

p = (^,0,0,/x) (B6) 

correspondingly. Then the subgroups of P(l,3) which keep (1B4I) . ( 1B5I) or ( 1B6I) invariant are 
the Euclid group E(3), the Poincare group P(l,2) in (l+2)-dimensional space or the extended 
Galilei group G(l,2) respectively. The basis elements of the corresponding Lie algebras are 

(-Pi) Pi, P3, J12, J23, J31), 

(B7) 

(-Plj P2, Pi, J\2, Joi, J02) 

or 

(-Pi) P2, P3, J12, Gi = Joi + ^31) G 2 = J02 + J32, D3 — D2 + ^03) (B8) 

where P a , J^ v and D 2 are generators given in ([9]) and (IBlj) . In addition, in all cases (1B4[) (1B6|) 
there are additional symmetries 

p ° = ik (B9 > 

and 

D\ = dpiiv . (BIO) 



The list of symmetries (1B7|) - (1B8j) had been found in paper 17| starting with conservation 
laws for equations ([3]) and ([5]). Our approach is more straightforward and guarantees finding of 
all Lie symmetries including nonvariational ones. An example of such nonvariational symmetry 
which cannot be found with approach used in [ijj is the dilatation whose generator D\ is given 
by equation (IBlOp . 
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